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The Dzyaloshinski-Moriya interaction partially lifts the magnetic frustration of the spin- 1/2 oxide 
SrCu2(B03)2- It explains the fine structure of the excited triplet state and its unusual magnetic 
field dependence, as observed in previous ESR and new neutron inelastic scattering experiments. 
We claim that it is mainly responsible for the dispersion. We propose also a new mechanism for 
the observed ESR transitions forbidden by standard selection rules, that relies on an instantaneous 
Dzyaloshinski-Moriya interaction induced by spin-phonon couplings. 



PACS numbers: 

Strontium Copper Borate (SrCu2(B03)2) is a new ex- 
ample of a magnetic oxide with a spin gap |lj], with a 
ground state well described as simply a product of mag- 
netic dimers in two dimensions on the bonds giving the 
strongest magnetic exchange M. The weaker exchanges 
are frustrated by the geometry and, as shown by Shas- 
try and Sutherland || , the ground state of the isotropic 
Hamiltonian is independent of the value of the weaker ex- 
change, up to a critical value. The excitations, however, 
are not purely local and cannot be explicitly given. Re- 
cent experiments by ESR Q and neutron inelastic scat- 
tering presented here show how in fact there are spin 
anisotropies needed for an accurate description of the dy- 
namics. We shall show the corrections to the ground state 
are needed that, while small, will be necessary to under- 
stand many physical properties. For example at finite 
external magnetic field SrCu2(B03)2 appears to exhibit 
a number of finite magnetization plateaux jj],|5| , and the 
anisotropies will determine the observability of plateaux 
in different field directions. Furthermore SrCu2(B03)2 is 
believed to be close in parameter space to a quantum crit- 
ical point whose nature is somewhat controversial, and 
while the anisotropies are small they may be essential to 
its nature. 

For spin | the leading anisotropic terms are of form 
Dzyaloshinski-Moriya || and exchange anisotropy. The 
former is particularly relevant, since it may not be frus- 
trated even if the isotropic exchange is. While a small 
Dzyaloshinski-Moriya interaction should not destroy a 
gap generated by larger isotropic interactions it modi- 
fies the pure locality of the ground state correlations and 
delocalizes the first triplet excitation. This is because it 
appears in lower order in perturbation theory than the 
frustrated isotropic interactions. In this paper we pre- 
dict the Dzyaloshinski-Moriya interactions that should be 
expected in SrCu2 (60*3)2 from the structure, and show 
that they do indeed explain new features of the excita- 
tions observed with ESR and neutron inelastic scattering 
experiments. 



Miyahara and Ueda [Q have introduced the frustrated 
Shastry-Sutherland model 



H 
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for SrCu2(B03)2, with S — 1/2 and where nn stands for 
nearest neighbor spins and nnn for next nearest neigh- 
bors. The lattice is shown in fig. |]. J = 85K and 
J' = 54A are antiferromagnetic interactions estimated 
from the susceptibility and the gap The spectrum 
of spin excitations has several interesting features [Q , in 
particular the existence of singlet bound states ||. The 
figure ||a summarizes the gaps to the first excited states, 
calculated by exact numerical diagonalization of the fi- 
nite size system (20 spins). In addition to the triplet 
state (solid line) calculated in , the energies of the two 
lowest singlet states (dashed lines) are given as a func- 
tion of J' j J. These roughly agree with recent calculation 
P): we find a ratio J'/J ~ 0.62, comparing |n| the cal- 
culated ratio of the energy of the Raman-active singlet 
(which is not the lowest singlet ||) to the triplet excita- 
tion with the same ratio determined by experiment . 
We find that the energy of the lowest singlet crosses the 
ground state at (J'/J) c ~ 0.68 (while the triplet state re- 
mains gapped) , in agreement with the transition recently 
predicted JO], but no evidence for the reported S = 1 in- 
stability y. Below this value, the ground state of (|l|) is 
simply a product of localized singlets. 

Anisotropic behavior of the first triplet energy, accord- 
ing to the direction of an external magnetic field, first ap- 
peared in ESR data Q , and cannot be explained by the 
fully isotropic model. We show that the Dzyaloshinski- 
Moriya coupling, which occurs in low-symmetry crys- 
tals, explains these features. In SrCu2 (60-3)2, an almost 
perfect center of inversion at the middle of the dimcr 
bonds forbids the Dzyaloshinski-Moriya interactions be- 
tween the two spins of a dimer. Each dimer is, however, 
separated from the neighboring dimer by a BO3 triangle 
for which there is no center of inversion at the middle of 
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the bond, allowing such an interaction between the spins 
of the nearest neighbor dimers. As the copper (ab) plane 
is approximately a mirror plane for the crystal structure 
the main components of D must be perpendicular to the 
plane thus lying along the c axis. We will neglect the 
other components, expected to be smaller. Furthermore, 
using mirror planes perpendicular to the copper plane 
and passing through the dimers (m\ and rri2 in the fig. 
|l|), we find an alternation of the D vector from bond to 
bond. Finally, the mapping of one dimcr onto the next 
one fixes the whole pattern of D vectors (fig. |l|). We 
therefore predict an anisotropic term to the Hamiltonian: 



(2) 



where Si and Sj are spin operators which belong to next 
nearest neighbors, the sign ± depends on the bond (see 
fig. [l]). e c is the unitary vector in the c direction. The 
magnitude can be estimated from Moriya's argument to 
be D ~ £±2±J' where Ag c ~ 0.28 has been measured 
by ESR @f and J' ~ UK §. This gives roughly D ~ 
0.5meV. 
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FIG. 1. Spin model including the Dzyaloshinski-Moriya 
interactions whose vectors are perpendicular to the plane. 
The arrows show the order of the spins in the expression 
D.(Si x Sj). At right, the unit cell of dimers A and B with 
the classical groundstate with J' = 0, and definition of the 
reciprocal space. 

The Dzyaloshinski-Moriya interaction is not frustrated 
because of the alternation of its vector; classically, in the 
absence of the exchange interaction J', the spins of the 
neighboring dimers would be perpendicular. Quantum- 
mechanically, the product of singlets on each dimer is 
no longer the exact ground state. To first-order the 
ground state has corrections proportional to two triplets 

T Ha^(b.c) s « itf 1 ^ 1 - t b X ' t t 1 )- These corrections are, 
however, small and we checked by exact numerical diag- 
onalization that the position of the critical point sepa- 
rating the two non-magnetic ground states is not appre- 
ciably changed. 

In order to calculate the effect of the Dzyaloshinski- 
Moriya interaction on the first triplet excitation, we start 



from small J' and from the exact ground state for D = 0. 
At J' = D = 0, the first excited state is a localized 
triplet which is separated from the ground state by an 
energy gap J. Because of the frustration, a finite disper- 
sion should appear only at the sixth order in J' / J [||. 
In contrast, the Dzyaloshinski-Moriya interaction is not 
frustrated and the degeneracy of the localized triplets is 
lifted to first order in D. We consider the two inequiv- 
alent triplet states for the dimers A and B (fig. 1). To 
linear order in J' and D, we can write the hamiltonian 
for the sector of total spin S z — +1 (resp. —1) as a 2 x 2 
matrix projected on the basis vectors: (i) All dimers are 
in singlets bar one, which is in the triplet state, spin +1 
(resp. — 1). This triplet is on the sublattice of dimers A. 
(ii) The same with the triplet on the sublattice of dimers 
B. One obtains: 
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where /(q) = cos(q a a/2) cos(<7(,a/2). The dispersion of 
the two modes ± (each is twice degenerate with S z = ±1) 
is therefore proportional to D: uj^(S z = +1) = u>^ (S z — 
— 1) = J ± 2D cos(q a a/2) cos(<7(,a/2). The Dzyaloshinski- 
Moriya interaction has no effect on the S z = compo- 
nent of the triplet, so that its energy remains equal to J 
[uJq(S z = 0) = J] (fig. |b). 
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FIG. 2. (a) Gaps to the first excited states at q = for 
D = (solid line triplets and dashed line singlets), calcu- 
lated by exact numerical diagonalization. (b) The effect of the 
Dzyaloshinski-Moriya interaction on the first triplet for small 
J', (c) For finite J', the exact diagonalization gives the ener- 
gies of two reciprocal points (the dots). (d)Renormalization 
of the splitting 8 and the bandwidths Wi and W%- 

In particular, at q = 0, we have two upper (resp. 
lower) modes S z = ±1 with lo(S z = ±1) = J + 2D 
[<->j(S z = ±1) = J — 2D]. A magnetic field parallel to 
z || D || c splits these modes in four branches. This is in 
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agreement with ESR when the magnetic field is parallel 
to the c-axis (fig. ||a). In such an experiment, however, 
an S z = state can not be seen as the external magnetic 
field is tuned to adjust the energy of the state to the fre- 
quency of the propagating wave. We have, however, a 
clear prediction for the energy of the mode S z — 0: it 
remains exactly at the middle of the four u>(S z = 0) = J. 
Such a crucial test can be made by neutron inelastic scat- 
tering experiments: we come back to this point next. 

A transverse magnetic field leads to the diagonaliza- 
tion of a 6 x 6 matrix which reduces to two equal 3x3 
matrices in the basis of the zero-field eigenstates. There- 
fore, the energies remain twice degenerate in a trans- 
verse magnetic field and are given by: = J ± 
^4L> 2 / 2 (q) + (. 9± /i S iJ_L) 2 , w° = J. This form also fits 
the ESR results in transverse magnetic field (fig. ||b), 
apart a small splitting in the high field regime which may 
be accounted by the differences in the g tensors from the 
dimer A and B. 




H (T) H (T) 

FIG. 3. The magnetic field dependence of the triplet ener- 
gies: (a) H || c, (b) H || a ; ESR data from j|, the dashed 
line comes from [fi3j. The theory of Dzyaloshinski-Moriya in- 
teraction explains these results (solid lines). 

These expressions are first-order in J'/ J (although cor- 
rections like DJ'/J 2 may appear), which is not accurate 
for SrCu 2 (B0 3 ) 2 (J 1 / J - 0.62). We performed exact nu- 
merical diagonalization to see how these results survive 
in the strong coupling limit. We know for instance that 
the gap between the ground state and the first triplet 
state is renormalized to order (J'/ J) 2 In addition, 
we find, for D = 0, a finite bandwidth W for the triplet 
dispersion. This bandwidth is smaller than that found 
previously from a perturbative treatment p2[ . We in- 
terpret this difference as breakdown of the perturbative 
series close to the quantum critical point. When D ^ 0, 
we calculated again the splitting 5 between lower and up- 
per modes at q = 0. The result S(J'/J — > 0) = AD is 
renormalized by finite J' (fig. ||d). Taking the relevant 
parameters, we find 6(J'/J = 0.62) ~ 2.01?. On the 
other hand, the dispersion is changed from that calcu- 
lated in the limit of small J': the energy of the S z = 
mode acquires a small dispersion W(D = 0) (fig. ||c,d), 
though smaller than the S z — ±1 dispersions. 



In order to test the prediction of the energy of the 
S z = mode, we compare to neutron inelastic scatter- 
ing measurements in the presence of a magnetic field, 
up to = 12 T. They have been realized on the 
three-axes spectrometer CRG/CEA-IN12 installed on a 
cold neutron guide at the ILL. The final wave vector 
was kept fixed at kf = 1.55 A -1 . Vertical curved py- 
rolytic graphite PG(002) monochromator and horizon- 
tally curved PG(002) analyser were used and a cooled 
beryllium filter was placed after the sample to eliminate 
higher-order contamination. The collimations were 40'- 
open-open. A single crystal of about 15 x 6 x 6 mm 3 
was installed in a cryomagnet with the c axis aligned 
along the vertical field. In the following, we refer to the 
(a*b*) reciprocal plane defined in (7) (fig. 0). In recipro- 
cal lattice units, the wave vector components are defined 
as Q a ,b = 9a,&a/(27r) where a is the lattice parameter. 
Examples of energy scans obtained at Q = (1,0,0) for 
H = and H = 6T are shown in fig. [I| In fig. ||a, we 
observe a broad single peak at E ~ 3meV, corresponding 
to the lowest triplet excitation branch (the signal above 
~ 4.4meV corresponds to a higher triplet jjj as well as 
the dashed line in fig. |^b). For H — 0, it is of note that 
the 3meV mode is appreciably broader than the energy 
resolution (~ 0.2meV at full- width half-maximum). In 
fig. the application of H shows a Zeeman splitting of 
this mode in five distinct lines. The five peaks become 
resolution limited. Their energies are shown in fig. 0a. 
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FIG. 4. Scans in the neutron experiment at Q = (1,0,0): 
(a) the broad signal comes from three superimposed modes 
(dashed lines), (b) Zeeman splitting of the modes. 

The energies of the S z = ±1 excitations at Q = 
(1,0,0) (solid dots) agree with the ESR field-dependent 
branches (Q = (0,0,0)), confirming the translational 
symmetry of the model. In addition, we find the en- 
ergy of the S z = mode of the triplet (solid squares) 
which can not be seen by ESR. It lies exactly at the 
middle of the four other branches, as predicted above. 
Furthermore, the observed intensities of the S z = ±1 
peaks are all equal and roughly 1/4 of that of the S z = 
peak (fig. [|b). We calculate the spin correlation function 
transverse to the scattering wave vector which gives the 
intensities. Because of the purely imaginary off-diagonal 
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elements in the matrix (||), the upper (S z — ±1) and 
lower (S z = ±1) have the same intensity at zero field. 
Both are split in longitudinal field into two modes of half 
the zero-field intensity. Similarly, for the S z = sec- 
tor, there are, in fact, two degenerate modes at zero field 
which remain degenerate at finite field with the original 
weight. Each weight is twice that of S z = +1 since only 
spin correlation functions transverse to the momentum 
transfer contribute to the intensities. The intensity fac- 
tor 1/4 then follows. The agreement with experiment 
suggests that the ratio of intensities is a good test of a 
mechanism based on couplings of Dzyaloshinski-Moriya 
symmetry, even in the regime of large J' / J. 

Since the spectra are well fitted for both directions, we 
can extract the coupling D from the zero-field splitting 
8 = 0.352 ± 0.008meV. Using the numerical result 8 ~ 
2.0D (fig. 2d), this gives D ~ 0.18meV, consistent with 
the approximate value from Moriya's formula, as given 
earlier. We obtain the gap (which is no longer simply J 
but the function shown in fig. |^a) and the ^-factors in 
both directions: g\\ = 2.24 ± 0.09 and g ± = 2.00 ± 0.09. 

We then discuss the origin of the finite intensities in 
ESR B and infra-red absorption (TJ] . No magnetic tran- 
sition can be induced from a purely singlet ground state. 
Spin-orbit couplings generate anisotropies that may per- 
mit transitions from the ground state to the original 
triplet excited states. First, note that Dzyaloshinski- 
Moriya interactions are not sufficient to explain the ob- 
served transitions. Indeed, the magnetic operator Hm = 
h J2i $t 1S °dd under the symmetry transformation: the 
reflection using the mirror plane mi and the rotation of 
7r around the dimer bond. The ground state and the 
triplet state are, however, even under this transforma- 
tion. The matrix element therefore vanishes, leading 
to zero intensity. Second, the exchange anisotropy only 
(both in J and J') leaves the ground state as the prod- 
uct of singlets, leading again to zero intensity. With both 
Dzyaloshinski-Moriya and exchange anisotropy, the tran- 
sition can exist but should be very small, of the order of 
A 6 , where A is the spin-orbit coupling. Note also that an 
exchange anisotropy of order D 2 /J always occurs with a 
Dzyaloshinski-Moriya coupling H , but unlike the unfrus- 
trated case [jlJJ, here it does not eliminate the splitting 
which is linear in D. Third, we have mentioned that the 
anisotropy of the g tensors may be responsible for the 
very small splitting of the modes in high field. A stag- 
gered field between dimers (due to perpendicular orienta- 
tions of oxygen ions surrounding copper ions) is written: 
H M = hAg (£ ieA £+ ~ E lG s St)- This term is odd 
in the transformation introduced above and transitions 
are thus forbidden for the same reasons. An intra-dimer 
staggered field should be very small giving an intensity 
we estimate as A 2 a 2 where a - a few degrees - is the 
tilting angle of the crystal structure |lq] . 

Consider now electric dipole transitions between mag- 
netic states. We restrict here to a phonon-assisted mech- 



anism. A purely electronic mechanism may in fact also 
apply Treating the spin-phonon interactions per- 

turbatively leads to an effective matrix element between 
pure magnetic states | (/ | He \ i)\ 2 , where the effective 
operator He has now a correction of spin-orbit origin 
: H e = E„„7Si.Sj + 7 ? .(S ; x Sj). The first term 
comes from a transition via a virtual phonon as an in- 
termediate state. The second term, which comes from 
the spin-orbit coupling should be | n |~ A7, with a di- 
rection that depends on details of the phonons. This 
mechanism gives an intensity proportional to rf, suggest- 
ing that the observed transitions may be electric dipole. 
Further polarized experiments can test this: if the in- 
tensity of the transition changes in rotating the crystal 
around the magnetic (resp. electric) field of the wave, the 
transition would be demonstrated explicitly to be electric 
dipole (resp. magnetic). 

We have shown that the static Dzyaloshinski-Moriya 
interaction explains the magnetic field behavior of the 
first triplet excitation (fig]|) and gives a dispersion 
(fig.||), because it is not frustrated. We extracted its 
value from the zero-field splitting, taking into account a 
finite J': D = 0.18meV. Since the direction of the main 
anisotropy is given by D || c, we expect clear magneti- 
zation plateaux only when the magnetic field is applied 
along the c-axis. This is consistent with current results 
[j5| . The static interaction does not, however, explain the 
finite intensity for the ESR transitions. We presented an 
alternative mechanism using a dynamical Dzyaloshinski- 
Moriya interaction. 

We thank Prof. H.Nojiri for providing us with his ESR 
data and for interesting discussions. 
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